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Abstract 

We clarify the hnks between a recently developped long wavelength itera- 
tion scheme of Einstein's equations, the Belinski Khalatnikov Lifchitz (BKL) 
general solution near a singularity and the antinewtonian scheme of Tomita's. 
We determine the regimes when the long wavelength or antinewtonian scheme 
is directly applicable and show how it can otherwise be implemented to yield 
the BKL oscillatory approach to a spacetime singularity. When directly ap- 
plicable we obtain the generic solution of the scheme at first iteration (third 

1 



order in the gradients) for matter a perfect fluid. Specializing to spherical 

symmetry for simplicity and to clarify gauge issues, we then show how the 

metric behaves near a singularity when gradient effects are taken into account. 
PACS numbers: 98.80.Cq, 04.50.+h 
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I. INTRODUCTION 

Despite the fact that the universe is clearly inhomogeneous on galactic scale and the 
possiblity, raised by some inflationary scenarios (see e.g.) [Q, that its geometry be "chaotic" 
on scales larger than the Hubble radius, most cosmological models are still the homogeneous 
and isotropic model of Friedman, Robertson and Walker (FRW). Many convincing reasons, 
physical or philosophical, can be given to that state of affair, but there is also a purely 
technical one : very few inhomogeneous solutions of cosmological interest are known (see, 

e.g. @). 

Various approximate solutions however have been given in the past. A simple one is the 
"quasi-isotropic" solution of Lifchitz and Khalatnikov (see e.g. [jlO|) the spatial sections of 



which (in a synchronous reference frame) are just uniformly stretched in the course of time 
[ds"^ = ~dt^ + a'^{t)hijdx'^dx^, where the arbitrary "seed" metric hij{x^) depends on space 
only). This metric is exact and reduces to the standard FRW metric if the spatial sections 
are maximally symmetric, and is a good approximation to an exact solution of Einstein's 
equations if, as we shall recall below, all spatial derivatives remain small, that is if all "point 
to point" interactions are neglected. 

A more general approximate solution when all gradients are neglected is the "antinew- 
tonian" solution of Tomita's [Q, which, as we shall recall, depends on as many arbitrary 
functions as a generic solution of Einstein's equations. 

Finally the "general oscillatory solution" studied by Belinski, Lifchitz and Khalatnikov 
is the most elaborate approximate description of a generic solution of Einstein's equations 
near the Big-Bang. 

There was recently a renewed interest in these approximation solutions first of all because 
the observations of the COBE satellite urged a fresh view on the old problem of structure 
formation (see e.g. |0) but also because a new line of attack of Einstein's equations was 
pursued. Indeed, in a series of papers , Salopek Stewart and collaborators developped a 
"long-wavelength" iteration scheme not of Einstein's equations but of the Hamilton- Jacobi 



equation for General Relativity. Their method, which consists at lowest order in neglecting 
all spatial gradients, leads back in most instances to the quasi-isotropic solution mentionned 
above, and yielded, for dust at least, the solution up to and including the third iteration 



(that is accurate to order 6 in the gradients as will become clear below |T2[)- (The case of 
a more general perfect fluid is more awkward to handle in this Hamiltonian formalism.) 

In another series of papers the present authors (together with Comer, Goldwirth, Tomita 
and Parry) [|l|, |^, |^ iterated in the same way the Einstein equations themselves. They 
noted that the zeroth order quasi-isotropic solution, although not generic, is an attractor at 
late cosmological times of the generic solution of Tomita. Concentrating then on this quasi- 
isotropic limit of the zeroth order solution they obtained the solution up to and including the 
second iteration (5th order in the gradients) for matter being a perfect fluid with constant 
adiabatic index or a scalar field. In the particular case of dust their result is identical to 
that of PI so that the link between the two methods could be clearly made. 

The motivation for going beyond the zeroth order is the wish to describe inhomogeneities 
within the Hubble radius. Indeed the approximation at the root of these long wavelength 
iteration schemes is the following. Take a synchronous reference frame where the line element 
reads: 

ds"^ = -dt"^ + -fij{x^,t)dx'dx^ {hh= 1,2,3). 

At each point define a local scale factor a and a Hubble time H^^ by 

a^ = (det7ijy/^ , H = a/ a, 

where a = da/dt. The Hubble time is the characteristic proper time on which the metric 
evolves. The characteristic comoving length on which it varies is denoted L : di'jjk ~ L^^^ij. 
The long wavelength approximation is the assumption that the characteristic scale of spatial 
variation is much bigger than the Hubble radius, that is: 

At lowest order then the long wavelength approximation is not suited to describe e.g. the 
formation of structure within the present Hubble radius. One can hope however that the 



iteration scheme pushed at a sufficiently high order can give results valid within the Hubble 
radius. Some numerical investigation of this question has been undertaken by Deruelle and 
Goldwirth (see also |p) but further work is nevertheless required to assess the convergence 
properties of the approximation scheme. 

Another motivation to go beyond lowest order, which is the one for this paper, is to 
study, within a long wavelength approximation scheme, the behaviour of a generic solution 
of Einstein's equations near a space-time singularity and make the link between that scheme 
and the BKL general solution referred to above. 

The point of view is therefore very different from the one adopted in our previous papers 
since, instead of the late time quasi-isotropic solution, we consider here the solution near 
a singularity, thus going backward in time. In this case the quasi-isotropic behaviour is no 
longer valid and the approximate solution becomes more complicated since it is no longer 
possible to separate the time dependence and the spatial dependence into a scale factor and 
a "seed" metric respectively. The generic (i.e. without assuming quasi-isotropy) ffist order 
solution was given by Tomita in the case of dust and radiation (it can also be found in |]^). 
Here we consider the more general case of a barytropic perfect fluid with an equation of state 
of the form p/t = F — 1, where F is a constant. Although an explicit solution for the first 
order solution cannot be given for F different from F = 1 (dust) and F = 2 (stiff matter), 
an explicit limit near the singularity can be given in all cases (sections 2 and 3). 

Once the first order solution has been given, we analyze the validity of the approximation 
near the singularity (section 4). To do this we examine the time evolution of the terms which 
were neglected at first order. We find that they could not always be ignored and we give a 
condition of validity for the approximation scheme. In the cases when this condition is not 
fulfilled, we are able to make the link with the work of Belinski Kalatnikov and Lifschitz. 
We believe that the way we recover the oscillatory behaviour of the metric, which does 
not introduce intermediate Bianchi IX geometries, is more straightforward than the original 
approach of BKL, and will allow in particular an easier analysis of the genericity of the 
"spindle" singularities found by Bruni et al [jTll] . 



We then give the generic third order solution (section 5). For the sake of simphcity we 
then apply in detail the approximation scheme to the case of spherical symmetry (section 
6). In particular we show that, in the case of dust, the third order solution corresponds 
to an expansion of the Tolman-Bondi solution in time. Finally, in section 7, we give our 
conclusions and comment on the usefulness of the long wavelentgh approximation. 

II. THE LONG WAVELENGTH ITERATION SCHEME 

In this section we first rewrite Einstein's equations in a way convenient for our purposes 
and then describe the iteration procedure. 

We place ourselves in a synchronous reference frame where the line element takes the 
form: 

ds^ = -dt^ + -fij{t,x^)dx'dx^ {i,j = 1,2,3) (2.1) 

(Coordinates transformations involving four functions of space can still be performed without 
spoiling the synchronicity of the reference frame; see e.g. |1^.) Matter is taken to be a perfect 
fiuid with pressure p, energy density e, unit four velocity u^ (/i = 0, 1, 2, 3) and stress energy 
tensor: 

T^w = (e + p)Uf,u,. + P9i,u (2.2) 

with the further restriction that p/e = F — 1 where the index F is supposed to be constant, 
positive and less than to 2 (the limiting cases F = and F = 2 correspond respectively to 
a cosmological constant and a "stiff" fluid whose speed of sound equals the speed of light; 
F = 1 is dust, F = 4/3 radiation; fluids with < F < 2/3 violate the strong energy condition 
and can be called "inflationary".) 

Einstein's equations are : W^ = Sj^ with Sj^ = xi^lf — ^6^Tp and x = SttG, G being 
Newton's constant. In a synchronous reference frame the components of the Ricci tensor i?(^ 
are (see e.g. [PH]]): 
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where Kij = 'jij is the extrinsic curvature (a dot denotes the derivative with respect to time 
t, a semicolon the covariant derivative with respect to '~fij); all indices are raised with the 
inverse metric 7*-'; k, = k] and i?* is the Ricci tensor associated with the metric 7jj. 
Now, one can always decompose /t* = Y^ikj into a trace and traceless part: 

4 = 2H5i + A'^la"" (2.4) 

where the "anisotropy matrix" A* is traceless (A* = 0) and where we have introduced a 
local "scale factor" a = (det7ij)^/®, so that the local "Hubble parameter" H is H = a/ a. 
When a{t,x'') and AUt,x^) are known the metric '-fij is obtained by integrating the 6 linear 
equations: 

7i, = 2H-f,j + 47fe*/a^- (2.5) 

(The matrix A* is therefore such that Af'-f^j = A'j'jki.) 

Let us then rewrite Einstein's equations as equations for a and Ay The traceless part of 
Kj = Sj gives: 

i} = 20=^(5] - i?}) (2.6) 

where Rj = i?* — |5}-R and Sj = xer(-u*-Uj — ^S^u'^Uk) are the traceless parts of i?* and Sj. 
The trace of /?* = S^ together with the (q) equation R^ = Sq give: 

2H + 3TH^ + ^r(2 - r) = 1(2 - 31)^ - ixe(4 - Sr)!^^ (2.7) 



Xe(l + r.z'^Wfc) = 31/2 - M! + i^ (2.8) 

where \\A\\^ = AiA{. 

Finally the (°) equation i?° = 5*° reads: 

2d,H-U^) =XeTu,^l + u''Uk. (2.9) 



Equations (|2.6| - |2.9|) are strictly equivalent to Einstein's equations but are written in a 
form suitable for the implementation of the long wavelength iteration scheme. We shall also 
use the following consequence of Einstein's equations (obtained by differentiating ( p.8| ) and 
using i^MM))- 






euitt yl + u^Uk 



(2.10) 



The long wavelength approximation consists in neglecting all the terms quadratic in the 
gradients, that is in the spatial derivatives, in Einstein's equations. Now, from Eq ( |2.9| ) or 
( [^.1U| ) the 3-velocity Ui is at least first order in the gradients so that the right-hand side of 
Eq ( |2.6D is at least second order. At first order then it can be set equal to zero so that Eq 
( [^.6D trivially gives that the anisotropy matrix A* does not depend on time: A* ^ ^^^ A^Ax^). 
Then Eq ( |2.7|) , the right hand side of which can be ignored at first order, is an equation 
which, when integrated, gives the scale factor *^^^a(t,a;^). The anisotropy matrix and the 
scale factor being known, Eq (|2.5| ) yields the first order metric *^^^7jj(t, x'^). Finally Eq (|2.8| - 
|2.9|) where at first order R and u^Uk can be ignored give the energy density ^^'e{t^x^) of 
the fiuid as well as its 3-velocity ^'''>Ui{t,x^) (the three velocity can equivalently be obtained 
from ( p.lO|) ). This first order solution is given by Eardley, Liang and Sachs [0, and Tomita 
p. It is reviewed in the next Section. 

III. THE GENERIC FIRST ORDER SOLUTION 



In this section we give the general solution (see JTSl and 0) of the truncated Einstein 



equations ( p.6[ - |2.9|) in which all terms of order greater than one in the spatial derivatives 
are neglected. 

A. The anisotropy matrix 

As already mentioned in section 2, Eq ( |2.6| ) at lowest order reads A* = and gives that 
the anisotropy matrix depends on space alone: 
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A] = (1)4 (x*^) with (^M^ = 0. (3.1) 



B. The scale factor 



As for the Eq { \i.7\ j for the scale factor a it reduces to: 



H + 3TH'^+^-^ ^(2-r) = 0, (3.2) 



a first integral of which is readily obtained: 



where a{x'') is an integration "constant" and where we have set /3^ = (3||*-^''A|p/8) . (We 



chose a > which will correspond to spacetimes emerging from a singularity. The collapsing 
situation d < is the time reversal of the solution presented here. As for a it will have to 
be positive or zero: see Eq ( p.l4| ) below.) Eq ( |3.3| ) can be explicitely integrated when the 
anisotropy matrix vanishes (/5^ = 0) or when matter is dust (F = 1), a radiation fluid 
(r = 4/3), or a stiff fluid (F = 2) (the particular case F = is treated at the end of the 
section); we give the expression for a in the case F = 4/3 for completeness only : 

r = 1 : (i)a=^ = uia^'u + (3^) (3.4) 



o /Q4 

r = 4/3: u = — r^ (xVa;2 + 1 - In (Va;^ + 1 + a;)) with x = 2d^^^^ a) / (3^ (3.5) 



r = 2 : ^^^a^ = c?u with a^ = ^13^ + 453 (3,g) 

where u = t — to{x'^) with to{x^) an integration "constant". For a general < F < 2 (and 
/3^ 7^ 0) an approximate solution for small u is: 

0<F<2 «a3 = /5\[l + £n2-r + o(n^-2r^] with d=^^r'^. (3.7) 



C. The metric 

The anisotropy matrix and the scale factor being known, Eq ( |2.5| ) then gives the metric 
^^^'jij{t,x''). Indeed, assuming that ^'^^A'^ is diagonahzable (this does not spoil the genericity 
of the solution), let us denote by r(a)(x''') its three eigenvalues and by Cq three associated 
independent eigenvectors with components el^{x^) that we shall normalize to unity. The 
triad ef forms a basis of the tangent space. Defining the cotriad e° by e°;el = 6^, the 
components of any tensor on this new basis is obtained by contracting its components in 
the coordinate basis with the triad or cotriad. In particular we define: rjab = 'jijel^el and 
Kab = f^ij^a^i- Since the triad is time independent the relation Hab = rjab holds and Eq ( p.5|) 
becomes: 'jab = f(^a)lab/o? with 7afc = rjab/o?'- Moreover the fact that A^-^jk = A^'jik implies 
that T]ab is diagonal (see |jl| for a demonstration). Hence the solution for the metric: 

^'hij = e'^ef'^Vab with ^^^ab = 5aba^e^v{r{a)jdta-^^ (3.8) 

where we recall that r(a) are the eigenvalues of the anisotropy matrix and e\ three associated 
normalized eigenvectors (the tracelessness of A* implies that Sr(a) = 0; we also have Sr? -^ = 
|p)yl||2 = 8/?73). Explicit integration of Eq (^ with a given by Eq (^ ^^ gives: 



r = 1 (i)r/,, = 5abC^a)l3^''^^'>'^"u"^+'-^^->/^\ua^ + l3yri~n',)IP' (3.9) 

r = 2 ^'\ab = 5ab-a'C^a)U^/-'+^(^^'-^' (3.10) 

< r < 2 ^'\ab = 5abP'/'C^a)u'/'^''^^^/^' 



' (") \ „,2-r I /n f„A-2r 



l + c~|^-....:^)^^-^ + 0(« 



3 /32(2 



(3.11) 



where C(a)(a;'^) are three integration "constants". We have that det7ij = a^ — 
det^(e")det?7afe = det^(e")C(i)C(2)C*(3)a^, and therefore det^(e")C(i)C(2)C(3) = 1. We note 
that when expanding the metric ( p.9| ) in small u, one finds the metric ( p.ll[ ) for F = 1. 



We also note that at leading order the metric ( |3.11| ), and therefore, because of the pre 



vious remark, also the metric (|3.g| ), is Kasner like. Indeed, setting p(^a) = 1/3 + r(a)/2/5^, 
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it reads ^^^rjab oc n^^^"' with Sp(a) = ^pfa) — 1- However it is important to note that the 
metric ( p.lOj ) for F = 2 is NOT Kasner-hke, since the sum of the pJ^^ (with the definition 
P(a) = 1/3 + r(a)/2a^) is less than 1 but one has still Sp(a) = 1. 

Finally, we should mention that the previous calculations are only valid for /?^ ^ 0. If 
/3^ = then the anisotropy matrix vanishes, and the integration of eq (|2.5|) is obvious. The 
metric is quasi-isotropic and reads 

^'hi, = a\t)h,,ix'), (3.12) 

where hij is an arbitrary "seed" metric that depends only on space and the scale factor must 
be taken from the paragraph b. depending on which type of matter one considers. This 
particular case of a quasi-isotropic metric was studied in details in our previous paper [Q] . 

D. The energy density 



As for the energy density it is given by Eq ( p.8[ ) which reduces to: 

Xe = 3i72-||(^)A||V8a^ (3.13) 

that is, using ( p. 3D : 

Xe=ld'a-'^. (3.14) 

(The positivity of e implies that a has to be positive.) For ^^^a given by Eq ( p.4| , pTB| - p.7D , 
(^) yields: 

^-' ^"^-^ 3.(.+V/a3) ^'-''^ 

0<F<2 x^^)e=^(3-'^u-^[l-5ru'-^ + 0{u^-'^)], (3.17) 

and one notices that the surface u = -^ t = to^x'^) is a singular surface of infinite density. 
(Note also that if, in the general case < F < 2, eq ( |3.14| ) yields ( |3.17| ), eq ( |3.13| ) only gives 
the leading part in xe.) 
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E. The three velocity 



As for the three velocity at first order it follows from ( p. 91) . However Eq ( 2.1(]| ), which at 
first order in the gradients reads 

~3^{3-3r)p^. = (1 _ r) jdta^diia^a-^^) + Q, (3.18) 

gives us the behaviour of Ui without having to resort to the full expression for the metric. In 
the case of dust it tells us for example that Ui is a function of space alone, and for < F < 2 
it gives: 

(^)t.. = d^to + a«^-^[l + (r - l)ln''^] + ^¥^, (3.19) 

with Ci = Cid~^P'^^^~^yT. To determine the three "constants" Ci{x^) as functions of the 
constants C(a) appearing in the metric, the more complete eq (|2.9|) must be used; their 
explicit expressions in the case of spherical symmetry will be given in section 7. 

F. Genericity 

Let us now examine the genericity of the metric thus obtained. It depends on the follow- 
ing 12 arbitrary functions: ^^^A^x'') (8 functions) (or, equivalently: r(a)(x^) (2 functions) 
and e^(x^) (6 functions)); a{x'') and to(x^); and the two functions C(a)(x'^). 

Now 4 of these 12 functions can in principle be fixed by choosing a particular synchronous 
reference frame (see section 7 for an explicit implementation of such a gauge fixing in the 
case of spherical symmetry). One sees in particular that the reference frame can be chosen 
in such a way that the surface of infinite density be t = 0, that is one can choose to to be 
zero. Indeed in an infinitesimal change of coordinates: i = t + T and x* = x* + X* with 
T = T(x*) to preserve synchroniticity, the three velocity transforms as Ui = -i^u^ + ^Uj ~ 
—diTuo + {6j — diX^)uj ~ Ui — diT if \ui\ << 1. Therefore the 3- velocity can be set equal to 
zero by an appropriate choice of coordinates if it depends on space only, which is the case 
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for dust. In all other cases this freedom of gauge can be used to set dito = as can be seen 
from Eq ( ]3J9| ). 



The metric ^^^'jij therefore depends on 12 — 4 = 8 physically distinct arbitrary functions 
of space corresponding to the 4 degrees of freedom of gravity (the two gravitons) and the 
four degrees of freedom of a fluid (e and m*). It is therefore generic. 

G. Late time limit 

The scale factor being an increasing function of time, Eq ( |3.3|) tells us that when a is 
large the anisotropy /3^ becomes negligible (unless F = 2) and the scale factor tends to 
its Friedmann- Robert son- Walker value a oc u'^^^^ ~ t^/^'" (tg can be neglected for large t). 
Since then J dta~^ oc t('"^2)/r _^ g^ -g^ (|3.8|) tells us that ^^^'jij tends to a "quasi-isotropic" 
metric : ^^^'jij — > f^/^^hijiyX^) where hij{x^) is a "seed" metric depending on three physically 
distinct arbitrary functions. Five physical degrees of freedom are therefore diluted away: 
the traceless part of the intrinsic curvature and the epoch of the Big-Bang (in the particular 
case F = 2, the metric does not become quasi isotropic at late times and only the epoch 
of the Big-Bang is lost). The quasi-isotropic scheme developped within an Hamilton- Jacobi 
framework by Salopek, Stewart and collaborators |§] and along the lines presented here 
in DH, which consists in iterating Einstein's equations starting from the restricted "seed" 
'"^^lij = t^^^^hij(x''), is therefore justified far away from a spacetime singularity. On the 
other hand, near a singularity, the full first order metric must be taken as starting point. 

H. The case of vacuum 

The Einstein equations for vacuum can be derived from the general equations for a 



perfect fluid by imposing in ( p.6[ - 12.9| ) e = and Sj = 0. Note that in this case all the terms 
in ( |2.7] ) proportional to F cancel because of ( p.8| ) with e = 0. 

In the vacuum case the constraint (^.81), that is ( |3.14| ), gives a = and Eq ( |3.3| ) gives 
(i)a3 = p'^u so that the metric (U) reads: T]ab = SabP^^^C(^a)u'^P(-^ with p(„) = i + r(„)/2/52. It 
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is a Kasner-like metric since Ep(a) = Sp? x = 1. (This incidently shows that matter becomes 
neghgible near u = 0: see Eq ( p.ll| ).) The metric depends on 11 functions but ( p.9| ) gives 
three additional constraints. Therefore the solution depends on 8 functions that is 4 physical 
degrees of freedom, those of the gravitons, as it should in vacuum. 

I. The case of a cosmological constant 

In the particular case of a cosmological constant (F = 0), the integration of ( p.3|) gives: 
/3V : (^^a^ = ^= sinh(V3An) , p^ = : ^^^a^ = exp{V3Au) (3.20) 



with A = Aa^/3 and the metric ( |3.8|) reads: 

^4\ 1/3 

3A 



/?V : ^''>Vat = Sa, (^] '" q.)(sinh y3A74t.)2/=^+^V)//5^(cosh ^^Auf'-^-^^^'^' 



VK 



P' = 0: (^, = 63 %.(x^). (3.21) 

(In the case /3^ = the integration constant to can be absorbed in the seed metric hij.) Eq 
becomes a definition of the energy density: 



X 



We = A = 'Vonst." (3.22) 



and ( 2.10 ) says that A is a true constant, independent of space, and hence is not a true 
degree of freedom. The metric therefore depends on 11 functions, 3 of which disappear 
when the constraint (|2.9| ) is imposed. In that case then, as in vacuum, the solution depends 



on 4 physical degrees of freedom, as it should. 

IV. CONDITIONS OF VALIDITY 

The purpose of this section is to establish in which situations the approximation scheme 
developped in the two previous sections is valid and which remedy to give in the cases where 
it is not. To check the validity of the approximation scheme, we simply compare the third 
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order terms arising from the first iteration which were ignored up to now, with the first order 
terms that we have just calculated. Far from a spacetime singularity when the first order 
metric reduces to its quasi-isotropic component this was already done in ^ and |1|] with the 
conclusion that the next orders tend to zero as time increases if matter violates the strong 
energy condition, i.e. if the fluid is "inflationary". On the other hand, near a singularity 
where the anisotropy matrix cannot be ignored, the next order, as we shall see below, blows 
up generically as one approaches the singularity, whatever equation of state matter satisfies, 
and we shall recover the BKL (Belinski- Khalatnikov-Lifshitz) oscillatory behaviour for the 
metric. 

Let us begin with the most tiresome part: the computation of the Ricci tensor built 
from the first order three-dimensional metric. We shall here compute the Ricci tensor for a 
general metric of the form 

1^3 = e-e'ji^ab, (4.1) 

where the triad e° depends only on the spatial coordinates whereas the metric rjab depends 
both on spatial coordinates and time. The situation is therefore more complicated than if 
rjab were only time dependent as is the case in the BKL analysis. However we proceed along 
similar lines. 



Let us first introduce the Ricci rotation coefficients (see 1JT4[); defined by 



7a6c = e(a)i;fce(b)e(^), (4.2) 



and their commutator. 



Kbc = 7a6c — 7ac6, (4.3) 



which have the property 



7abc = -Ibac + dcVab, (4.4) 



thus enabling us to express the Ricci rotation coefficients in terms of their commutators: 
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labc = \ \Khc + \ca - Kab + dcVab + d^rjac - daTjbc] ■ (4.5) 

It is important to rewrite the commutators Xabc in the form 

Aa6c = Vadl^^c + ^cVab - dbTjac, (4.6) 

where 

Ac = e^e^ (d,ei - d^ef) , (4.7) 

are time independent. In order to make the link with the terminology of BKL it is worth 
noticing that these coefficients can be rewritten in the form 

It is then straightforward to compute the components of the three dimensional Riemann 
tensor in the nonholonomic basis. One finds 

Rabcd = VecVfbdal/^ - VecVfadblf^ + Idebl^ca ' Ideal^cb + ldcel\a - Idcel^b' (4-9) 

Therefore the components of the Ricci tensor are 

Rab = Vebdair ' V^'VebVfadcl/^ + l\,l\, - ^bel^,. (4.10) 

In terms of the functions fiabc, the Ricci rotation coefficients can be expressed as 

labc = \ \^abc + ^bca " ^cab + dc^lab + daTlhc " Qb^la^^ ■ (4.11) 

SO that the explicit expression of the Ricci tensor in terms of the triad and of the metric riab 
is given by 

Rab = dafl\^ - ^dcfl\a ' hVebV"^dcfi\d + lVeaV"^dcfi'db 

+ 5/^ ec (/^ ba + t^ba ~ t^a b) ~ 1 it^ ae ~^ t^ae ~ t^e a) (/^ be + t^bc ~ t^c b) 

+ ht^\cV^'^ (ddVab - daVbd - dbTlad) - IfJ'^aV^^ deVcf + I (/if,/ - f^ab) V^^ C^dfVec " OeVfc) 

+ \lJi'"^a C^dcVbd + ddljbc - dbTfcd) + IjJia'^dd'nbc + ll^'^^'b (deVaf " ^a^e/) " ^jlb^'dfriea 

_^l^cd^e/ (^d^riaeddTlbf - dcVbedfrjda) + Iv""^ {daddVbc - daOb^ldc - OcOdVab + dbOcVad) ■ (4.12) 
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In addition to the curvature terms we have also ignored, in the first order approximation, 
the terms quadratic in the three-velocity. The approximate three- velocity that follows from 
the first order metric according to ( p.9| ) is given in the new basis by 

XeVua ^ -\ (d,K\ + k\^\, - nil',, - 8^^) , (4.13) 

where we have used the relation 

4,e^e5e^ = d^K^, + ^7^ - <l\c- (4.14) 

We have now to analyze the time dependence of these two expressions giving the Ricci 
tensor and the three- velocity and find the dominant terms, i.e. that with the smallest power 
in u since we are heading towards the singularity. To do that we shall assume that 1 . the new 
metric is diagonal, i.e. rjab = Va^ab (this amounts to supposing that the anisotropy matrix is 
diagonalizable) and 2. the spatial derivative of any component of the metric has the same 
time behaviour than the component itself (this means in particular that we assume that to is 
independent of space - which, as we saw, can be the case in an appropriate reference frame 
-, and that we also ignore the logarithmic corrections that arise from the spatial derivative of 
the exponents). Inspection of the above expression then shows that all the terms in Rat can 
be classified in one of the following categories, as far as their time behaviour is concerned: 
cste (or logarithmic), rja/rjc, rjh/rjc, rjc/rj^ and rjaVb/VcVd, where a and b are fixed but c and d 
range from 1 to 3, whereas the three-velocity, because simpler can be given explicitly here: 



Xe 



T^'^Ua - -h 



^ Vbh Vah ^ r. Vb 

l-b^a—1 ab 1 ab- ^b^aOa — 

Vb Va \Vb, 



(4.15) 



with 



l\b = fJ'\b + \ (dbVac + dcTjab " daTjbc) ■ (4.16) 

A. General case 

Now, as was shown in section 3, the generic behaviour of the first order metric near the 
singularity is of Kasner type, when < F < 2. Let us label the coordinates in such a way 
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that pi < P2 < Ps- We know that -1/3 < pi < 0, < p2 < 2/3 and 2/3 < ps < 1 (see 0). 
In Rii, the dominant term then is 

l(/23)%i/^2r/3 (4.17) 

(there is no term of the form rjl/rjl because of the antisymmetry of /i"^^, in the two last 
indices). In i?22; the dominant term is 

-l{^^\^)%l/m (4.18) 

whereas the dominant term in i?33 is 

-\{^^\^)\lh2. (4.19) 



As for the dominant contribution in the crossed terms, /2i2, -R23 and -R31, it is more compli- 
cated since they are several terms involved. Therefore we quote only the time dependence, 

R12 ~ ^i/^/s, ^23 ~ est, ^31 ~ rii/ri2. (4.20) 

The dominant term in the scalar three curvature R is 

-f(/23)'^- (4.21) 

The three curvature thus behaves as a power law, 

R ^ „2(pi-p2-p3) _ y4pi^-2^ ^4 22) 

As we see the dominant terms always come from the cross-product of the ^abc so that the 
terms with the spatial derivatives of the metric rjab do not play a role near the singularity. 
One can therefore expect that we will recover the results obtained by BKL who started their 
analysis on a Bianchi IX model where the metric rjab is only time dependent. 

The time behaviour of the Ricci tensor being now known, let us see if we were allowed 
to neglect it. An analysis of the Einstein equations ( |2.6| ), ( |2.7| ) and ( p^.Sj ), shows first that 



u Rl and Ru must be convergent for the approximation to be valid. This is NOT the case 



if the dominant terms are those hsted above because pi is negative. Indeed u^R is divergent 
as well as u'^Rf and u'^R\. 

As for the term containing the three- velocity in the equation (|2.7| ), it goes like 



u^ieVu^Ua) ~ u^-^""'. (4.23) 

Therefore this term is convergent only if F > 2p^. Moreover, since — 1/3 < pi < and 
2/3 < p3 < 1, one can see that when this term diverges it can be either more (F < 4j9i + 2^3) 
or less divergent than the curvature term. 

The conclusion therefore is that the long wavelength approximation scheme breaks down 
in the general case when approaching the singularity. 

B. Case /i^23 = 
In the particular case where 

li\, = 0, (4.24) 

all the dominant contributions listed above vanish and the validity of the scheme must be 
reconsidered. In that case the dominant time dependences are the following: 

-Rii ~ Vi/V3, R22 ~ ^2/^3, -R33 ~ est, (4.25) 

R12 ~ ^2/%, ^23 ~ est, ^31 ~ est. (4.26) 

Therefore 

^ ~ 1/773 ~ n-^P^ (4.27) 

Knowing that 2/3 < P3 < 1, one can conclude that Ru^ is then always convergent. In a 
similar manner one sees that all the quantities u^Rl are convergent. 

Let us now look at the term quadratic in the velocity. The time dependence of this term 
remains the same as in the general case. Therefore the approximation of equation ( p.7|) is 
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valid only (except the case F = 4/3) if F > 2^3 or in the case where M3 = if F > 2p2 or 
U2 = 0, these conditions being rather restrictive. 



This is however not yet the end of the story. The analysis of eq (|2.8| ) imposes two further 
conditions: 

1. u°'Ua must be small with respect to 1, which implies that F > 1 + Pa for all a unless 
Ua vanishes. 

2. R must be negligible with respect to the energy density e: this is due to the fact, 
already mentionned, that the two first terms on the right hand side compensate at 
leading order. The condition is therefore, in view of equations (|3.15[ - |3rT7D , that u^R 
must be convergent. Therefore one must have F > 2^3. 

Note that the condition F > 1 +^3 implies all the other conditions, but this is very restrictive 
in general since ^3 is limited from above only by 1. 

To summarize, we find that the first order solution given in section 3 is a good approx- 
imation to a solution of Einstein's equations near a spacetime singularity, if the conditions 
yU^23 = 0, F > 2p3 and T > 1 + pa (for all a = 1, 2, 3 unless Ua = 0), are satisfied. Diagonal 
anisotropy matrices form an important subclass of matrices satisfying the first condition and 
we shall see that in the context of spherical symmetry, studied in detail in the last section, 
the other conditions are also fulfilled for F = 1 or for F > 4/3. Now imposing ([4.24|) renders 
the first order solution non generic as it then depends on 7 instead of 8 physically distinct 
arbitrary functions. However a qualitative analysis of what happens when /i^23 7^ can be 
given, which follows closely the work of BKL. 

C. BKL oscillatory behaviour 

We now consider the Einstein equations ( |2.6[ - ^.7D where we do NOT neglect the curvature 
terms any longer. The time derivative of equation (P^), after use of ( |2.7| ) and ( p^.8| ) and 
expressed in the new basis, gives 
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a-' 



(a3<) • = (2 - r)xeS-, + 2xeru'^u, - 2RI (4.28) 



As shown in section 3, the energy density e goes hke u~^ . We now assume that the curvature 
term evolving hke u^^^~^ is dominant over the velocity term, evolving like vF~'^~'^^^ (at worse). 
Assume moreover that the metric remains diagonal in the evolution and that at some time it 
is Kasner-like. In these conditions we are able to recover the behaviour discovered by BKL, 
initially in the case of Bianchi IX and extended later to inhomogenous situations. 
Following BKL, let us introduce a new time defined by 

r = lnM, (4.29) 

and let us write the metric in the form 

[r^ab] = Diag[e^'^,e^^,e^'^]. (4.30) 



Then keeping in ( |4.28| ) only the dominant contributions from the curvature term, we get the 
three following equations governing the coefficients of the diagonal metric, 

a^a = -i(/23)'e^^ ^lP = \{^^\,)"e'-, dl^ = l{l^\,)\'\ (4.31) 

When the metric is Kasner-like, one has 

dra = pi, dr/3 = P2, dri = P3, (4.32) 

The equation for a is similar to the equation for a particle with coordinate a moving in an 
exponential potential. Initially the particle moves with a constant velocity dr<y = pi- After 
refiexion on the potential wall, the particle will move with the velocity dr<y = —pi- The two 
other equations then give the two other final velocities pi: dr(3 = p2 + 2pi and drj = p3 + 2pi. 
Therefore the initial Kasner-like metric evolves into another Kasner-like metric due to the 
infiuence of the curvature terms, given by 

r]ab ~ Diag[u^+^pi , u i+^pi , u i+^pi ] (4.33) 

We thus recover from our general analysis the oscillatory behaviour between Kasner-like 
metrics, behaviour which was studied in details by BKL (see e.g. |T( 
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Let us conclude this section with the particular cases of the vacuum, a cosmological 
constant and a stiff fluid. In the case of vacuum, one has still a Kasner-like metric as was 
shown in the subsection 3.h. Therefore the above analysis applies without modification. The 
same conclusion arises from the cosmological constant case with anisotropy since the metric 
( see subsection 3.i), when u ^ has the same Kasner-like behaviour as the metric (^.11[). 



However the case F = 2 gives a qualitatively different result. Indeed, there is more freedom 
for the coefficients p(^a) and it is then possible that pi, the smallest of the three powers, be 
positive, in which case all the terms u'^Rj converge. Note that the energy equation is valid 
also only in the case where pi > 0. BLK showed that in the case of stiff matter (or scalar 
field) the oscillating behaviour (if there exists one at the beginning, i.e. if pi is negative) will 
end after a few oscillations when one goes backwards in time. This means that, sufficiently 
near the singularity, the approximation scheme works. 

V. THE GENERIC THIRD ORDER SOLUTION 



Let us first rewrite the Einstein equations (|2.6| - 1279|) in the new basis: 

At = 2a%S^-Rt) (5.1) 



2H + 3TH^ + ^^(2 - r) = 1(2 - 3T)R -ix^(4-3r)r^.<'u. (5.2) 



Xe(l + Tu'^Ua) = 3H^ - ML + 1^ (5.3) 



where \\A\\^ = ^t^i^ and 

XeTua = - J-^ (d^< + ^\l\, - <l\, - da^) ■ (5.4) 

In the previous sections we have only considered the first order approximation of the 
Einstein equations. We can now include third order corrections to the first order quantities: 

a = (i)a + (=^)a + . . . , A^ = (^M^^ + (3)^^ + . . . . (5.5) 
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These third order corrections follow from the approximate Einstein equations in which the 
terms that were ignored previously are now taken into account but are computed with the 
first order solution. Since the third order solution is supposed to be small with respect to the 
first order solution, one can linearize the Einstein equations and all the equations giving the 
third order terms will be linear ordinary differential equations. One sees that, in principle, 
one can repeat this procedure at any order and build iteratively the metric and the other 
quantities. 

It is convenient to define the following third order quantities: 

as^—^ (5.6) 

and 

where a boldfaced letter stands for a matrix. The third order correction to the extrinsic 
curvature, '•^•'Ka, can be decomposed, following ( |2.4|) , into 

^'^Kl = 2^'^HS'^ + {yst (5.8) 

where 

^^^H = ds. (5.9) 

The expansion of Einstein equations then gives 

(3)^^ = 2 /dt Wa3(^fe«-i?«), (5.10) 

where the term under the integral is built from the first order solution and not from the 
exact solution as is the case in the exact Einstein equation ( |2.6| ) (Note also that the constant 
of integration that could arise from the above integral is supposed to be already included in 
the first order term *^^M||J). As for the equations for a and e they become: 

2d3 + GTHas + ^^^^ (^) .ys = ^i? - ^^^^^x^'hu^u.. (5.11) 
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The third order correction to the metric can be computed by expanding formula ( p.4|) 
(in the new basis). One finds 

^'^«^ = ^""^acV"' - VacV^'^'^VdeV'"'- (5-12) 

Using the fact that the matrix [rjab] is diagonal, 

^'\ab = Diag[ri,], (5.13) 

one finds 

where there is no summation on the indices. This equation can be integrated into 

%ab=Va[dt'^^'^Kl (5.15) 

J Va 



and therefore 



(3)^. =^ frl^!Iif0(.3)Tr;^b 



%b = Va dt^[2^'^H6'^ + {ys):}. (5.16) 

One can then obtain the third order energy and velocity by inserting the metric ^^^r]ab + ^^^'i]ab 
in equations ( ^.3|) and ( |5.4| ) (taking for the quadratic term UaU^ the first order solution ^-^^Ua). 
For instance the third order energy density is given by 

;^(3)g = e^iJag - -W ( — ] -ya +hR-xr^'^eWu-Wua. (5.17) 

VI. THE EXAMPLE OF SPHERICAL SYMMETRY 

A. The equations 

The line element of a spherically symmetric spacetime can be written in a suitable 
coordinate system [t, r, 6, (p) as: 

ds'^ = -dt^ + -irrdr"^ + leedn'^ with d9? = dO'^ + sin^ ed(f)'^. (6.1) 
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An infinitesimal change of coordinates t = t + T,r = r + R preserves the synchronicity of the 
reference frame if T = T{r) and R = T'/'jrr (where a prime denotes a derivative with respect 
to r). It involves two arbitrary functions of space: T{r) and the integration "constant" in 
the equation for R (see e.g. |p!0| ). 

The extrinsic curvature k* = Y'^'ijk is diagonal and therefore the matrix A^- in Eq (4) 
is too, so that we are spared from the triad formalism of sections 3-5. We shall give here 
the solution to third order in the gradients. This will illustrate the general discussion of 
the preceeding sections, will clarify the gauge issue and allow a comparison with the known 
exact solution of Tolman-Bondi. This section is intended to be self-contained. 

When the line element is ( |6.1|) so that the traceless anisotropy matrix Aj in k^ = 2H6j + 
Aj/a^ (where H = a/a) is diagonal with eigenvalues r(r) = —A, r^e) = r^^) = A/2, the 



Einstein's equations ( p.6[ - p^ for a perfect fiuid yield 



7r 



a^ exp — dt 



A 



a^exp 



dt' 



A 



(6.2) 



where A and a are given by: 



A = -2 dta\S-R) with S 



2xe 



Vu^'ur , R = Rl-\R 



(6.3) 



2H + 3ri/2 + A ("4 ) (2 - r) = 1(2 - 3r)i? - i(4 - 3r)^. 

16 \a-^ 



(6.4) 



As for the energy density and the radial velocity they are given by: 



Xe(l + rnV)=3if2-^(^) +\R 



(6.5) 



Xer-Uryl + Tu'Ur 



^ (li 



A 



4V^^)(ln7ee)'. 



(6.6) 



Simple counting gives that a generic metric solution of ( |6.2[ - |6.4| ) depends on 4 arbitrary 
functions of r. The total number of physical degrees of freedom however is 2 (there are 
no gravitons in spherically symmetric spacetimes and the fiuid is specified by its density 
and radial velocity). Two functions can therefore be eliminated by fixing the gauge that 
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is choosing a particular synchronous reference frame, in agreement with the remark below 
Eq ( |6.1| ). We can first give a geometrical meaning to the coordinate r by relating it to 
the surface of 2-spheres: this will fix a function in 'jee- To eliminate the remaining gauge 
freedom we note that in an infinitesimal change of coordinates that preserves synchronicity, 
the radial velocity transforms as: Ur = ^Uq + ^Ur — —T'uq + ( 1 — ^j m^ ~ m^ — T' if 
Ur « 1, so that an arbitrary function of space (T') can be substracted to u^- 

B. The case of dust (r = 1) 



At first order in the gradients Eq (|673|) gives A = A{r) and the solution of (|6.4|) is: 



^^'a^ = a^u{u + a) with u = t — to and a = |-^ (6.7) 



^1 

3~ 



a 



and where a{r) and to{r) are 2 integration "constants". The metric then follows from ( |6.2| 
— see eq (p.9|): 



where e = A/|A|. From now on we shall consider only e = 1 since the metric, in the 
case e = —1, tends toward a metric of the type Kasner with the coefficients {pi = l,p2 = 
0,ps = 0) which is nothing less than the fiat metric as can be shown with a suitable change 
of coordinates (see ||10|)- therefore there is no singularity for e = —1. Cr and Cg are 2 
integration constants. The first order metric depends, as anticipated, on 4 functions of r: 
a,to,Cr and Cg and 2 can be eliminated by fixing the gauge. To do that we first impose 
Cg = r'^, so that the radial velocity ( |6.6|) becomes: 



(Dm, = t'o + ^a (6.9) 

which, as we already knew from ( p.l9| ) depends on space only and can be set equal to zero 
(indeed when matter is dust and hence follows geodesies there exists a synchronous reference 
frame where the particles remain at rest) by choosing a = — 2tQr/3. 
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The gauge being thus completely fixed, the first order spherical symmetric line element 
for dust finally reads: 

Wds^ = -de + -^^dr^ + o'dQ'' (6.10) 

1 — hr^ 

with the definitions: q = rv?!"^ ^ u = t — to and C^ = (1 — hr'^)^^ (this form will be useful 
for a comparison with the exact solution of Tolman and Bondi). It depends on 2 arbitrary 
functions of r: h and to- The energy density follows from ( |6.5|) (see Eq (|3.15|) ) and the radial 
velocity is zero: 

X^'^e = ^ ^ , , , , ^'^^^ = 0. (6.11) 

Useful secondary quantities are: 

(i)a3 = ~aMu - 2t'r/3) ; (^^ ( ^) = , "^^"^ , , . (6.12) 

^ ° ' ' ' Uy 9u{u-2t'or/3) ^ ' 

To obtain the third order metric the easiest way is as follows: Writing a = ^^^a(l + 
as); A/a^ = ^^\A/a^) + y^, Eq (IJ) gives: 

^^Us^ + ^^Us^ = -de + ^, J l + -ir)dr^ + q\1 + 7e)d(]' (6.13) 

1 — hr^ 



with 7r = 2a3 — Jdtys, 'jg = 2a^ + i Jdiy^. Eq (|6^ ) for 03 and 1/3 is then transformed, by 
using the relation y^ = —3Hy-^ — 3'^^)(y4/a'^)d3 — 2(5 — R) which follows from (|6.3| ), into an 
equation for 7^ which eventually reads: 

7e + 37e [^^'^H + l^'H^/^')] -R+IR = 0. (6.14) 

With ^^^H = ^^^d/a and ^^\A/a^) given by (|6.12| ), and the relevant components of the Ricci 
tensor for the line element (|6.1CI|) being: 

- (hry 2hr'^ ~ 2hr'^ 2ihP)' 

R = ^—-, :—^ ; R=^^ + —, (6.15) 



Eq ( |6.14| ) becomes 7e + 2'jg/u + hu ^^^ = which is readily integrated into: 

^e = -^hu'/\ (6.16) 
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Turning then Eq ( |6.4]) into an equation for 7^. instead of 70 and using the relation 7^ = 7e — fz/s 
to ehminate y^ one gets: 



7r 



+ 37. l^^^H - l^^\A/a^)] +2R+IR- l^e^^\A/a^) = (6.17) 



the integration of which yields: 



g .,5/3 c 2/3 

^^ = -^(h + h'r) '^—— + -hrt' ^^— . (6.18) 



At third order then, the spherically symmetric element for dust is ( |6.13|) with 70 and 7^ 
given respectively by (|6.16|) and ( |6.18|) . 



Now the exact solution for dust with spherical symmetry is known. It is the Tolman- 
Bondi solution, the line element of which can be written as (see e.g. Landau Lifschitz [0): 



with 



for / > 0, 



ds^ = -dr + — ^TT^dr^ + p^dVL, (6.19) 

l + /(r) 



P= ^(coshr/-l), t - to(r) = ^^(sinh r/-r/) (6.20) 



p= ——{l-cosvi), t-to{r) = ———^{T]-smT]) (6.21) 



for / < 0, and 



P = Y^\t-to{r)f' (6.22) 



for / = 0. The line element is written in the comoving gauge where the three velocity 
vanishes and the miscalleneous functions appearing in the metric have an easy physical 
interpretation: any particle is labelled by the coordinate r, the same at any time; 47rp^(r, t) 
gives the area of the sphere containing this particle at time t; p{r, t) is the radial velocity of 
the particle and fi{r) corresponds to the mass inside the sphere containing the particle. 

We now consider the expansion of the Tolman-Bondi solution in the parameter u = 
t — to{r), supposed to be small. We find for p: 
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p 



ox 2/3 Q /On 2/3 



'fu'/'. 



(6.23) 



In the derivation of this expansion, we assume that u' is of the order of u. Going to the next 
to leading order in the expansion we thus find : 



7r 



1 /3\^/=^ {fxu 



9 V2 



2\2/3 







1 + / 



(in/iM^j 



1 + 



3/^2/3 /2\V' 



(in fiu^j 



ilnfu'/fi) 



(6.24) 



5/x2/3 V3; V"^^- J (ini^u^)' 

Identifying /(r) = —hr"^ and /i = Ar^/d, we recover the first and third orders ( |6.13 , 5.16 ) and 
( p.l8| ) given by the expansion scheme. 



C. The general case (0 < T < 2) 



The first order metric is obtained as before and, without going into details, it should be 



clear that it is given by ( ^.111 ), the index (a) being r and 6, with 



'(r) 



-A 



-4e/?V3;r(,) 



r((^) = A/ 2 (again we shall consider only e = +1 since e = — 1 does not describe a singularity). 
This metric depends on 4 functions of r: CrP^^^ = g{r), CqP^^^, to and c and 2 of those 
can be eliminated by particularizing the reference frame. As before we shall first choose 



CeP^^'^ = r^. Then, from (|3.19| ) and the discussion in section 3.f, we know that to can be 



chosen to be zero, and that fixes the gauge completely (note that when F = 1 it is not the 
comoving gauge chosen in the preceeding paragraph). The first order line element therefore 
with: 



IS 



(1) 



7r 



gt 



-2/3 



3 (2 - r) V 



(6.25) 



^'hee = rH'l^ 



1 + 



25(1 



j-2-r 



3 f2-r) 



+ o( 



4-2r 



(6.26) 



which depends on 2 arbitrary functions: gir) and cir). As for the density and radial velocity 



they are given by (|3.17|) and ( |6.6|) : 



;^(i)e=§(3-r)gt 



sj^-r 



i-art"-^ +c f 



i4-2r 



(6.27) 
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'"'- - ('"^)'r^' ^ r(3^;^ [' + ^(r - i)*-T . (6.28) 

Let us now examine the conditions of validity of the approximation by determining the 
behaviour of the third order terms. The analysis here is just a very particular case of the 
general discussion given in section 4. Indeed we are in the case where /i^23 = 0, U2 = U3 = 
and pi = —l/3,p2 = Ps = 2/3. Therefore we know that the curvature terms do not cause 
any trouble, that the condition UqM" < 1 implies F > 2/3, and finally that one must have 
r > 4/3 in order to insure the validity of equation ^]5| giving the energy density. The reason 
for which ( |6.27|) fails to give the energy density when F < 4/3 is that the leading order in 



( |6.27| ) compensate (the Kasnerian metric is a vacuum solution), and that the subdominant 
term is superseded, when F < 4/3, by third order terms coming from the first iteration. Let 
us recover these results directly by determining the behaviour of the third order correction. 
At leading order the dominant term in the r.h.s. of (|6.3| ) is ^ ~ —-^t~^/^; as for 
S it remains negligible {S oc t'""^/^). The first order scale factor being at leading order 
proportional to t^''^ we have that the third order correction to A, built out of the leading 
part of the first order solution is: 

which has to be compared with ^^\A/a^), built with the more accurate first order solution 
i ^lWM ), that is: 



(^) (^) = -^ [1 - ct'-' + O (t--)] . (6.30) 

We therefore see, in agreement with the general discussion of section 4, that indeed the first 
order solution is a good approximation to a generic solution of Einstein's equations near 
a spacetime singularity up to and including terms in t^~'" provided that 2/3 > 2 — F that 
is F > 4/3. When F < 4/3 the metric ( |6.25| - |6l26| ) is still good at leading order near the 
singularity but the energy density cannot any longer be given by (|6.27| ). 
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D. The case of stiff fluid (F = 2) 



The first order metric is again obtained by particularizing the results of section 3 to the 
case of spherical symmetry. The anisotropy matrix depends on r only: ^^^A = A(r); the 
scale factor is given by ( p.6|) : ^^^a = aw^/^, u = t — to; and the metric ( p.lO| ) becomes: 



(1) 



J TV \-^ f' Hi 



(1) 



7ee = L^eu 



^ 2/3+2a/3 



(6.31) 



(with a = SA/Aa^). It depends on 4 functions: Cr,Cg,to,a. We fix the gauge by choosing 
Cg = r^ and to = so that the generic first order line element, together with the energy 
density, given by (|3.16|) and the radial velocity derived from ( |OD are: 



Wds' = -de + t2/3 \c,.r^^i^ + rh'^'/^dQ' 



(6.32) 



X'^'e 



3t2 



(1) 



Ur 



1-a^ 



aa' In t + 



3a + ra' 



(6.33) 



(6.34) 



The solution depends on the 2 arbitrary functions: Cr and a, and the positivity of e imposes 
a2<l. 

To obtain the metric at third order in the gradients we must first evaluate the r.h.s. of 
equation (|6.3|), that is compute R and S by means of the metric ( |6.31| ). These are sums of 
terms in fr^^"^") [Const., Int, (Int)^] and in t~3(^+"), that we shall denote collectively by 
Q. Integrating (^ and (^ will then yield (3)^" = M + Qt^ and (3)a" = "(i)a(l + Qt^), 
so that the third order metric, given by (1) will be of the form: 



(l)+(3). 



Irr irr\^ ' ir) 



^''^^^'h9e = '^'hee{l + le) 



(6.35) 



where the time dependence of 7,. and 751 is Qt^. Hence we see that the iteration scheme is 
valid if Qt^ < 1. Since a^ < 1 this condition is satisfied near t = 0, in agreement with the 
general discussion of section 4. 
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The detailed calculation gives, for a > and at leading order near the singularity: 

9(1 — 2a) 2/r, „\ , , 

r^(7 — 2a){2 — ay 

(For a; < the dominant third order term near the singularity is in (lnt)^t3(^+"-'.) One can 
also check that the conditions of validity coming from the three-velocity are also satisfied 
since F = 2 and pa < 1- 



VII. CONCLUSIONS 

In this work, we have studied the early time behaviour of inhomogeneous spacetimes near 
the singularity. To do that we have used a long wavelength iteration scheme to approximate 
the Einstein equations. Our main concern was to test the validity of the approximation 
scheme, by comparing the terms we ignored with those we kept. The result of this investiga- 
tion is that one should be very cautious with the use of the long wavelength approximation 
if one wishes to get general results. Indeed our analysis shows that, in the general case, 
there are very severe restrictions on the range of validity of this scheme. The troubles arise 
from two origins: 

• The curvature terms: the curvature terms, which are ignored in the first step, blow 
up in general near the singularity. Going beyond the long wavelength approximation 
by keeping them from the beginning enables us to recover the oscillatory behaviour 
discovered by BKL. In the case where the curl of the vector field representing the 
axis of contraction (going forwards in time) is orthogonal to the vector field, then the 
curvature terms can be ignored and the approximation scheme is valid. 

• the velocity terms: the velocity terms may also blow up. 

They do not if the perfect fiuid is sufficiently "stiff" to compensate for the dilatation 
(F > 2pa) or if the component of the velocity along the dilatation (time going forwards) 
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direction vanishes. Even if these conditions are fulfilled it is not yet enough to get a 
valid energy equation, which demands that T > 1 + pa unless the component «„ 
vanishes. 

In view of these results it should be interesting to reconsider the study of BKL when the 
velocity terms are dominant over the curvature terms. In this case the role of matter should 
become important and one should not be able to restrict oneself to the case of vacuum. 

Finally we stress the fact that the problem of the velocity terms disappears in the case 
of a cosmological constant (F = 0), and in the case of irrotational dust where it is possible 
to choose a synchronous system of coordinates for which the three- velocity of dust is always 
zero. If we impose spherical symmetry the complete scheme works for F > 4/3 (including 
F = 2). However the scheme works weakly, i.e. without the energy relation, for F > 2/3. 
The case F = 1 is special since one can choose a coordinate system so that the scheme works. 
The scheme also works in general for the stiff case as soon as there is local expansion along 
all the spatial directions. 

All our conclusions should apply to a gravitational collapse, instead of a Big-Bang, by 
just reversing the time. 
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